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Abstract. We study the clustering properties of galaxy clusters expected to be observed by various forthcoming surveys both in 
the X-ray and sub-mm regimes by the thermal Sunyaev-Zel'dovich effect. Several different background cosmological models are 
assumed, including the concordance ACDM and various cosmologies with dynamical evolution of the dark energy. Particular 
attention is paid to models with a significant contribution of dark energy at early times which affects the process of structure 
formation. Past light cone and selection effects in cluster catalogs are carefully modeled by realistic scaling relations between 
cluster mass and observables and by properly taking into account the selection functions of the different instruments. The results 
show that early dark-energy models are expected to produce significantly lower values of effective bias and both spatial and 
angular correlation amplitudes with respect to the standard ACDM model. Among the cluster catalogs studied in this work, it 
turns out that those based on eRosita, Planck, and South Pole Telescope observations are the most promising for distinguishing 
between various dark-energy models. 



1. Introduction 

One of the main quests of contemporary astrophysics is the de- 
termination of the nature of dark energy, the mysterious com- 
ponent of the cosmic fluid responsible for the accelerated ex- 
pansion of the Universe. While the eviden ce of its presence 



has became compelling in the past decade ( Astier et al 



Riess et all 



Rubin et al 



2007 , iDunklev et alj 20081: iKomatsu et al. 



2006; 



2008 



2008tlKilbinger et all2008l) . its nature remains en- 



tirely unexplained. In particular, while virtually all the present 
observational evidence is in concordance with a cosmoiogicai- 
constant interpretation of the dark energy, its possible dynami- 
cal evolution is not well constrained. The detection of this time 
evolution would hint at dark energy being different from vac- 
uum energy, and would call for some m ore general explanation , 
such as minimally coupled scalar fields dRatra & Peeb les 1988; 
Wettericd[l988rjl iBrax & MartirfeOOOh . Since, in this standard 
generalization, dark energy is not supposed to clump on the 
scales of the largest cosmic structures, the only way its na- 
ture can be unveiled is by studying the expansion history of 
the Universe. The expansion rate of the Universe as a function 
of cosmic time in turn affects the process of structure forma- 
tion, and consequently the many observable properties of cos- 
mic structures that are accessible to observations. 

The most immediate effect of cosmology on cosmic struc- 
tures is on the number counts of objects, especially the most 
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massive and extreme ones such as galaxy clusters. A generic 
dynamical dark-energy model that postulates a dark-energy 
density increasing with redshift will necessarily imply an ear- 
lier structure formation than in cosmoiogicai-constant models 
(if the linear amplitude of density fluctuations at present is 
fixed), and thus a higher abundance of objects at any time. 

An alternative channel for the detection of possible effects 
of dynamical-dark energy is the study of clustering properties 
of cosmic structures. Models predicting a higher abundance of 
objects imply that high-mass clusters are less exceptional, and 
as a consequence the linear bias with respect to the underlying 
dark-matter density-field should be reduced. The linear bias, 
the structure abundance, and the linear correlation function of 
density fluctuations all affect the angular and spatial correla- 
tion functions that are observed in cluster catalogs, and all de- 
pend on the behavior of dark energy. We are therefore justi- 
fied in exploring the effect of different quintessence models on 
the clustering properties of massive galaxy clusters, and under- 
standing whether differences between them could be detected 
significantly in cluster catalogs produced by forthcoming ex- 
periments. This is the purpose of the present work. 

We focus on blind cluster surveys based on the X-ray emis- 
sion of the hot intra-cluster plasma and on the spectral distor- 
tion of the CMB radiation produced by the thermal Sunyaev- 
Zel'dovich ( Sunvaev & Zeldovichll972 , SZ henceforth) effect. 
The cosmological models that we address range from the con- 
cordance cosmoiogicai-constant ACDM cosmogony to early 
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dark-energy models, to intermediate models with a constant 
equation-of-state parameter for dark energy w x t - 1 or with 



Sadeh et all d2007b per 



a gentle evolution in w x with time 
formed a simple preliminary study of the two-point angular 
correlation function predicted to be observed by Planck in 
models with early quintessence. As will become evident from 
the results shown in this paper, our findings are qualitatively 
consistent with those, while a quantitative comparison is not 
directly possible because of the different catalog definitions. 
Observational results about the clustering properties of galaxy 
clust ers measured in opti c al and infrared ca talogs can be found 



m 



Brodwin et al. 
d2008h . 



(2007J), iPapovichl (120081) and Estrada et al 



This paper is structured as follows. In Sect. [2] a brief 
overview of the different cosmological models employed in this 
paper is presented, with a description of their main features and 
a summary of the various cosmological parameters. In Sect. [3] 
we review the formalism used to describe clustering of galaxy 
clusters in the past light cone of the observer. In Sect.Hl we de- 
tail the properties of the X-ray and SZ surveys analyzed in the 
present work and summarize the scaling laws used in linking 
the mass and redshift of objects to their observable properties 
in Sect. [5] In Sect. [6] we describe the properties of the clus- 
ter catalogs obtained therefrom. In Sect. [7] we present our re- 
sults on the spatial and the angular correlation functions, and in 
Sect. [8] we summarize our conc lusions. We shall use through- 
out the ISheth & Tormenl d2002l) prescription for the computa- 
tion of both the cluster mass function and the linear bias. 



2. Cosmological models 

We use seven different cosmological models. The first four 
of them are models with an early dark-energy (EDE hence- 
forth) component, labeled from EDE1 to EDE4. In early dark- 
energy cosmologies, the dark-energy contribution is assumed 
to be represented by a quintessence scalar field whose evolution 
tracks that of the dominant co m ponen t of the cosmic fluid at a 
given time (IWetterichll 1988allbl 11995). As a consequence, the 
density parameter for dark energy at very early times does not 
vanish as in more conventional models, but flattens to a finite 
value. To ensure dark-energy dominance at low redshift how- 
ever, an ad hoc mechanism for breaking the tracking behavior 
must be adopted, usually in the form of a non-standard kinetic 
term i n the quintessence Lagrangian (Hebeck er & Wetterichl 
200ll iDoran et al.l 12001b boran & Robbers! 120061) or a non- 
minimal coupl ing between qu intessence and neutrinos with 
evolving mass JWetterichll2007l) . 

An adequate parametrization of early dark-energy models 
consists of the dark-energy density parameter at present, O x ,o, 
the dark-energy equation-of-state parameter at present, w x ,o and 
a suitable average of the dark-energy density parameter at early 
times during the phase of linear structure formation, 



^x,sf = 



In 



i r° 

— Q. x (a)d\m 



(1) 



where a eq is the scale factor at matter-radiation equality. 
Observational constraints from Type la supernovae, large-scale 



Table 1. Parameter values for the seven cosmological models 
investigated in this work 



Model 




£2x,o 


"x,sf 


h 


w x ,o 


0"8 


EDE1 


0.330 


0.670 


0.040 


0.670 


-0.928 


0.820 


EDE2 


0.360 


0.640 


0.040 


0.620 


-0.997 


0.780 


EDE3 


0.284 


0.716 


0.033 


0.686 


-0.942 


0.715 


EDE4 


0.282 


0.718 


0.048 


0.684 


-0.935 


0.655 


K08 


0.279 


0.721 




0.701 


-1 


0.817 


w x = -0.8 


0.279 


0.721 




0.701 


-0.800 


0.817 


ACDM 


0.279 


0.721 




0.701 


-1 


0.817 



structure, an d CMB allow Q x s f to be o n the order of a few per- 
cent at most dDoran et alj|20 05. 2007) . Our models EDE1 and 
EDE2 were introduced and studied in Bartelmann et al. I (l2006l) 



(see also lFedeli & B artelma nnll2007ah. while EDE3 a nd EDE4 
were investigated bv lWaizmann & Bartelmann! ([2008) and have 
cosmological paramete rs more closely adapted to the latest 
WMA P data releases (IDoran & Robbers 2006; Doranetal 



2007). For a detailed analysis of how early dark-energy mod- 
els compare with other dark-energy models on observational 
grounds, especially with respect to type-la supernovae data 
sets, we refer the reader to lRubin et al. [ d2008h . 

Apart from the EDE models, we also investigate a cosmo- 
logical model with d ynamical dark-energy parametrised as in 
Komat su et al. (2008), which we briefly describe below. In this 
case, the dark-energy equation-of-state parameter is assumed 
to evolve with the scale factor as 



W x (fl) 



aw a 
a + a* 



a{\ — a)w\ 
a + a, 



a + a» 



(2) 



where z* = 1/a* - 1 is a transition red shift that we set to be 
Z, = 10 in what follows. In any case. iKomatsu et al.l (|2008) 
showed that the precise choice of the transition redshift is not 
extremely relevant to the inferred value of the parameters wq 
and w\. In the low-redshift limit, z <k Z*, Eq. (01 reduces to the 
more standard form 

w x (a) = w + (1 - a)wi (3) 

dChevallier & Polarskil20"oillLinderl2003l) . We shall assume in 
the following that wq = -1, while w\ = 0.5, which is almost 
the highest possible val ue inferred at 95.4% confidence level 
by WMAP-5 year data (IKomatsu et al.ll2008b . For brevity, we 
refer to this model by K08 in the remainder of the paper. 

In addition to those described above, we analyzed a model 
with constant w x = -0.8 and a standard ACDM cosmological 
model with parameters given by the latest WMAP-5 year data 
release in combination with Type-la supe rnovae and baryon 
acoustic oscillations ( IKomatsu et al.l 12008). The redshift evo- 
lution in the equation-of-state parameters for our seven dark- 
energy models is shown in Fig. [T] while the values of the main 
cosmological parameters are summarized in Tablefl] There, the 
Hubble constant is expressed as Ho = h 100 km s _1 Mpc -1 , 
and cr 8 represents the rms of primordial density fluctuations 
smoothed on a scale of 8 h comoving Mpc. 

Note the highly non-trivial behavior of w x (z) in the early 
dark-energy models, for which the equation-of-state parameter 
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1 2 3 

Log(l + z) 

Fig. 1. The redshift evolution of the dark energy equation of 
state parameter w x for the seven dark energy models adopted in 
this work, as labelled in the plot. 

reaches positive values for relatively low redshift, especially 
for EDE3 and EDE4. We also note the extremely low normal- 
ization of the power spectrum for linear density fluctuations 
(parametrized by <xg) in model EDE4, necessary to counter act 
the quite high dark-energy density at early times, which in turn 
determines a particularly low linear density contrast threshold 
for spherical collapse. As can be seen from Fig.Q] all our mod- 
els approach the cosmological-constant behavior at low red- 
shift, and they are all constructed to also be in agreement, 



beyond CMB data, with large-scale structure (Teg mark et al 
l2004allbl) and Type-la supernovae data. 

3. Clustering formalism 

We used the formalism developed by Matarrese et alj (19971) 



and f urther applied, among others, by iMoscardini et al 



( 1998b in studying high-redshift g alaxy clustering and by 



Moscardini et al. (2000, 2001, 2002) in describing the cluster- 



ing of galaxy clusters in the past light cone of an observer given 
a survey selection function. In this section, we briefly summa- 
rize this formalism and refer to the quoted papers for additional 
detail. 

The starting point of the formalism is the number of ob- 
jects of mass M per unit redshift around z given a back- 
ground cosmology, N(M,z) - 4-ng(z)n{ M,z), where n(M,z) 
is the standard differential mass functio n ( Press & Schechtei 
1974tlBond et al.lll99ll: Isheth & Tormenll2002l) and g(z) is the 



Jacobian determinant 
dr 



g(z) 



r 2 (z)^(z). 
dz 



(4) 



Evaluating the differential mass function for cosmological 
models with dynamical dark-energy, we set the linear-density 



threshold for collapse as computed in lBartelmann et aL ( 20061) . 
We shall further comment on this choice in Sect. [8] In Eq. (0), 
r(z) is the comoving radial distance to redshift z, and g(z) 
thus represents the comoving volume per unit redshift around 
Z. Equation (|4]i is valid only for a spatially flat cosmological 
model, which we assume henceforth. Under the same assump- 
tion, the comoving radial distance r(z) is 



r(z) 



c 



dz' 
E(z') 



(5) 



where E(z) = H(z)/Hq is the normalised Hubble parameter. 

We consider a galaxy-cluster catalog covering the mass 
range [Mi, M{\, where M\ and Mi in general depends on red- 
shift. Then, the true all-sky equivalent redshift distribution of 
objects in the catalog reads 



N(z)= I N(M,z)dM. 
Jm, 



(6) 



For the realistic situations explored below, the observed red- 
shift distribution equals Eq. © multiplied by the fractional sky 
coverage / s k y of the survey. In the real world, cluster catalogs 
are restricted by a threshold value of some observable (for in- 
stance the X-ray flux or the SZ decrement), hence all objects of 
mass higher than the mass M\ , corresponding to this limiting 
observable at a given redshift by some scaling relation, will be 
included in the catalog. This amounts to setting M2 = +°° in 
Eq. ©. 

One important ingredient for predicting the observed clus- 
tering properties of galaxy clusters is a relation between the 
density contrast of collapsed objects and that of the underlying 
matter distribution with the correct time evolution, i.e., the ef- 
fective bias. For the linear 'monochromatic' bias, we adopt the 
expression 



1 



b{M,z) = 1 + — 



Dl(z)S(M) 



Oc 



1 



l + [^6j(D + (z)^/S(M))] 



2/' 



(7) 



where D + (z) is the linear growth factor for density fluctuations, 
S (M) is the variance in the primordial density field smoothed 
on a scale corresponding to the mass M, and 6 C is the linearly 
extrapolated density contrast at collapse for a spherical density 
perturbation. We note that in general S c depends on redshift 
(although mildly for models with constant w x ), but for clarity 
we omitted this dependence in the previous equation. 

The standard iPress & Schechterl (1 19741) (see also 
Mo & Whitd ll996) relation is recovered by setting a — 1 and 
p - . The more precise relation p roposed bv lSheth & Tormen 



(1 1999b 1 see also lSheth et al.ll2001h can be obtained by instead 
setting a = 0.707 and p = 0.3. We shall follow this second 
option for consistency with the mass-function prescription and 
because it has been shown to provide closer agreement with 
the bias measured in ACDM numerical simulations. 



4 



C. Fedeli et al.: Clustering with dynamical-dark energy 



Given the above, the effective bias is defined as the linear 
bias weighted by the abundance of clusters in the catalog at 
hand, 



1 r A 



b(M,z)N(M,z)dM. 



(8) 



We define £(r, z\ , Zi) to be the two-point correlation func- 
tion of the underlying density distribution, for density peaks 
placed at the two different redshifts z\ and Z2- This is con- 
veniently computed by Fourier-transforming the non-linear 
power spect rum of the density fluctuations, described e.g., us- 
ing the fit o flPeacock & Dodds ( 199q) . More accurate prescrip- 
tions for evaluating th e non-linear matter power spectrum ex- 
ist (Smith et al. 2003), but their differences to that given by 
Peacock & Doddsl (1199a) are small and relevant only on scales 
< 1 Mpc h~ x , which are not pertinent here. The object correla- 
tion function can then be defined as 



&bjfazi»&) = b eS (zi)b eS (z2)^(r,ZuZ2) ■ 



(9) 



The problem given by the presence of a double redshift de- 
pendence in the correlation function of the underlying den- 
sity fluctuations is solved by considering a single, average 
redshift z defined as D + (z) = ^D + (z\)D + (z2) so that, effec- 
tively, %{r,Z\,Z2) = %(r,z)- The effect of redshift-space dis- 
tortions is also taken into account by mult iplying the c orre- 
latio n function £(r,zi,Z2.) with the factor (jKaiserl 119871 see 
also lZaroubi & Hoffmadll996UMatsubarall2000l) 1 + 2B(z)/3 + 
B 2 (z)/5, withjS(z) = f(z)/b eS (z) and 



f(z) = 



d\nD + (z) 
dMl+zY 



(10) 



In models with a cosmologi cal constant, th e function f(z) can 
be well approximated by ( Lahav et aDll991 ) 



f(z) 



Q 0.6 ( } ^ 



1 + 



(ID 



In more general models with a dynamical evolution in the dark- 
energy component, f(z) must be evaluated numerically. 

Taking all past light cone effects into account, the observed 
spatial correlation function is given by 



I r z »«t r z > 
= 42 I I 



'£inf -£irrf 

with normalisation 



N(zi) N(z 2 ) 
r{z\) r(zi) 



-'-in I 



r(z) 



■dz- 



€ob)(r,Zi,Z2)dzidz2, (12) 



(13) 



In the previous two equations, z; n f and z sup are the minimum and 
maximum redshifts, respectively spanned by the cluster catalog 
at hand. In realistic situations, zm - 0, while z sup depends on 
the sensitivity of the instrument used. 

Accordingly, the observed angular correlation function is 



Wobs(0) 



where 



•-'lint >JZmt 



N(z\)N(z2)^{T, zi , Z2)dzidzi, (14) 



f = Kzuzi, 0) = V^Czi) + r 2 (zi) - 2r(zi)r( Z2 ) cos(0) (15) 



and the normalisation in this case is the total number of objects 
included in the catalog at hand, 



B 



N(z)dz. 



(16) 



In the small angle approximation (see for instance iPeebles 

1980), the relation in Eq. ( fT4l simplifies to 



Wobs(#) 



1 r /V(-) r +M 

m TTTH ^ bi(r t ,,z,z)dudz, 

B J z . f dr(z)/dz J- m 



n mI dr(z)/dz 
with 

r, = r m {u,6,z) = Vw 2 + r 2 (z)9 2 



(17) 



(18) 



With the help of this formalism, it is possible to produce 
realistic theoretical expectations for the correlation properties 
of galaxy clusters in the past light cone, as a function of catalog 
properties and the cosmological model. 

4. Survey properties 

We consider five ongoing and planned cluster surveys here, two 
of which are X-ray based, while the remaining three are SZ 
based. 



4.1. X-ray catalogs 

Two forthcoming X-ray surveys are addressed in this work. 
The first has the properties of the eRosita wide survey, de- 
scribed in the mission definition document Q. It is designed 
to have a sky coverage of / s k y - 0.485 and a limiting flux of 
Fij m = 3.3 x 10~ 14 erg s~' cirT 2 in the energy band [0.5,2.0] 
keV. We note that this planned survey fulfills almost exactly 
th e requirements specifi ed in the dark-energy task force paper 
of lHaiman et al.1 d2005l) . where a survey of X-ray clusters opti- 
mal to constraining the evolution in the dark-energy equation- 
of-state is described. This proposed survey covers Q. 2 x 10 4 
square degrees (/ s k y 0.485) like the eRosita wide survey, and 
has a limiting flux in the [0.5, 2.0] keV energy band that is only 



slightly lower of, Fn m = 2.3 x 10 erg s cm . We adopt 

however the exact parameters of the eRosit a survey. 

The second is the XMM cluster survey (Sah len et al . 2008, 



XCS henceforth), a serendipitous search for X-ray clusters in 
the existing exposures of the XMM satellite archive. The sky 
coverage estimated on the basis of the already surveyed area, 
the pointings still to be analysed, and the expected mission life- 
time is Q. a 500 square degrees (/ s k y 0.012). As for the depth 
of the survey, since the XMM archive occupies a range of dif- 
ferent exposure times, and thus of limiting fluxes, a single lim- 
iting flux is inappropriate for describing the XCS catalog. If 
fluxes in the [0.1,2.4] keV energy band are considered, a sin- 
gle limiting flux of Fn m = 3.5 x 10~ 13 erg s~' crrr 2 might be 
used, although the redshift-dependent cut reported in Fig. 9 of 



Sahlen et al. 



Flim(z) 



10~ 13 erg s _1 cm 



(2008) is more appropriate, being fairly fit by 

= 2.8z~ 1/3 . (19) 



1 available at 

http : //www . mpe . mpg . de/pro j ects . html#erosi ta 
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We note that the XMM cluster survey has a brighter limit- 
ing flux at all cosmological redshifts and a much smaller cov- 
ered area than the eRosita survey, hence large differences are 
expected in the number of clusters, and the clustering signal- 
to-noise ratio detected in the former survey will be much lower 
than in the latter. 



4.2. SZ catalogs 

We consider three planned blind sub-mm surveys. The first is 
based on Sou th Pole Telescope (SPT) o bservations, and was 
considered bv lMaiumdar & Moh in an attempt to pre- 

dict possible constraints on cosmological parameters and dark 
energy. The proposed survey area amounts to Q 4 x 10 3 
square degrees (/ s k y - 0.097), for a limiting SZ flux density 
(see definition in the Sect. [5] below) of S Vo ,iim - 5 mJy at a 
frequency vo = 150 GHz. We note that a blind survey of SZ 
clusters with SPT has indeed already started with the first suc- 
cessful detections ( Staniszewski et al.ll2008l) . hence a compar- 
ison of our theoretical predictions with observations may be 
imminent. 

The second SZ cluster survey is based on the portion of 
sky that the Atacam a Cosmology Teles cope (ACT) will ob- 
serve. According to ISehgal et al.1 d2007l) . this will consist of 
two stripes covering 4 degrees in declination and 360 degrees 
in right ascension, for a total of Q. =s 3.6 x 10 3 square degrees 
(/sky — 0.087). The depth of the survey is defined in terms of 
the (frequency-independent) integrated Compton-y parameter 
over the solid angle covered by the virial sphere of each indi- 
vidual cluster, 7200 - According to the simulations performed in 



for z > 0.1 1, and by 



Sehgal et al. (2007), for a limiting integrated Compton param- 
eter of F200,iim - 10~ 3 arcmin 2 , a galaxy cluster-sample ~ 90% 
complete is produced, even with the inclusion of noise from 
radio and infrared point sources. 

The last survey that we consider will be carried out by the 
Planck satellite. This is the largest sub-millimeter survey of the 
sky that is currently being developed. Even though the beam of 
the satellite will be quite large, it is predicted to detect several 
thousands of clusters by thei r thermal SZ effect. As shown in 
Schafer & Bartelmann ( 2007b . the sky coverage for the detec- 
tion of galaxy clusters will be highly non-uniform, especially 
at low detection significance. On the other hand, we believe 
that a uniform sky coverage of Q 3 x 10 4 square degrees 
(/sky - 0.727) is realistic and sufficient for our purposes, and 
hence we adopt it here. As for the limiting SZ observable, the 
noise due to Planck's scanning path is highly structured on 
cluster-scales and shorter length. This means that assuming a 
simple flux-detection thresho ld is insufficient for our study. In 
Scha fer & Bartelmann ( 2007b . the minimum mass detected as 
a function of redshift is presented for a limiting integrated y- 
parameter of Jaxuim = 10~ 3 arcmin 2 , based on a numerical 
simulation in a ACDM world model. This minimum mass is 
virtually independent of the filtering scheme adopted, and is 
described reasonably well by 



log 



10 15 M R /i- 



-1.924 + 8.333z 



(21) 



if z < 0.11. This kind of fit may appear cumbersome, but we 
were unable to find a simpl er functional form tha t adequ ately 
reproduces the results of Schafer & Bartelmannl (2007), be- 
cause of the steep increase in the limiting mass at z > 0.1. 
Note that the two branc hes of the fit join smoothly at z = 0.1 1. 

The simulation by ISchafer & Bartelmannl d2007b was per- 
formed for a single cosm ological model, namely a WMAP-1 
cosmology (ISpergel et alj|2003b with Q m o = 0.3, Qa,o = 0.7, 
h = 0.7 and cr 8 = 0.9. To correct for the different cosmolo- 
gies used here, we proceed as follows. First, we convert the 
minimum mass from lSchafer & Bartelmannl d2007b into a min- 
imum integrated Compton y-parameter, according to the scal- 
ing relation described in Sect. |5]below and by considering the 
appropriate WMAP-1 cosmology. We then convert this min- 
imum integrated y-parameter back into a minimum mass us- 
ing the same scaling relation but using the various cosmologies 
adopted here. This procedure explains why clusters of the same 
mass produce different signals in different cosmologies, be- 
cause of their different formation histories and the differences 
in geometry of the universe between models. On the other hand, 
altering the cluster abundance may also change the amount of 
undetected objects, and thus the background noise. A proper 
treatment of this issue would require a far more detailed analy- 
sis and probably fully numerical simulations, which we decide 
was unnecessary for our purposes. 

5. Scaling relations 

To relate survey properties to the extent in mass and redshift 
space of the resulting cluster catalog, it is necessary to link the 
mass and redshift of an individual cluster to the relevant ob- 
servable, namely X-ray flux for X-ray surveys and SZ flux den- 
sity or integrated Compton-y parameter for sub-mm surveys. 
We do this by means of realistic scaling laws. We note that not 
all features of these scaling relations are well established, es- 
pecially concerning their redshift evolution. We described in 
the following what we propose to be the most suitable way to 
proceed given our aims. 

5.1. X-ray scaling relations 

First of all, we used the conversion bet ween the X-ray tern 



perature and the vi rial mass adopted in Fedeli & Bartelmanrj 
d2007bb (see also iBartelmann & White! l2003h . i.e.. a virial 
relation with normalisation based on the simulations of 



Mathiesen & Evrardl d2001l) 

M200 



kT(M 20 o,z) =4.88keV 



10 15 M C 



Hz) 



2/3 



(22) 



where the mass is measured in units of M @ . Additionally, a 
luminosity-temperature relation given by 



log 



Aftim(z) 

10 15 M o /r 



-1.200+ 1.469 arctan[(z - O.IO) 0440 ] (20) L(T) = 2.5 x 10 43 erg s" 



h- 



kT 



1.66keV 



2.331 



(23) 
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was used. This is based on observations by lAllen & Fa bian 
(11998b , and is assumed not to evolve with redshift according 



to the analyses of Mushotzky & Scharf ( 1997), Reichart et al 



to the analyse 

(1999) and lHashimoto et al.l (120021) . Combining these two rela 
tions, we obtain the mass-luminosity scaling law 



L(M 20 o,z) = 3.087 x 10 44 erg s^ 2 



M 



200 



10 15 M F 



■Hz) 



1.554 



(24) 



Choosing a reasonable value for the Hubble constant, h — 
0.7, Eq. (ED equals 



L(M 200 ,z) = 1.097 x 10 45 erg s" 1 



-E(z) 



, (25) 



where the mass 



is 



-l 



now expressed in M B h 
(2003) demonstrated that this mass 



B artel mann & W hite 
luminosity relation is a good fit to the X-ray c luster observa- 
tions compiled bv lReiprich & Bohringerl (12002b . 

A possible steepening of the luminosity-temperature rela- 
tion for low -mass clusters and groups of ga laxies has also bee n 
advocated dHelsdon & Ponmanl l2000albt IXue & Wul |2000|) . 



which would require replacing Eq. J23l with a broken power- 
law. However, it has been shown dOsmond & Ponmanl 12004 



Khosroshahi et al. I l2007h that the scaling relation for groups is 
consistent with that for clusters, although the scatter for groups 
is considerably larger. This could bias the estimate of the re- 
lation's slope; However since this has not been established 
definitively, we prefer to adhere to Eq. (f23t in the following. A 
steepening of the luminosity-temperature relation for groups of 
galaxies would have the consequence of including fewer low- 
mass objects in the various catalogs. These would thus con- 
tain a higher fraction of high-mass clusters whose clustering 
properties would enhance the differences between cosmolog- 
ical models. Our results may thus slightly underestimate the 
distinguishing power of the correlations. 

To convert the bolometric X-ray luminosity provided by 
the scaling relations to the luminosity in a given band re- 
quired to characterize the clust er catalogs of eRosita an d XCS, 
we adopt a Raymond-Smith dRavmond & Smithll 19771) plasma 
mode l implemented via the xspec software p ackage (|Ar naud 
19961) . with metal abundance Z = O.3Z dFukazawa et al.l 



1998; Schindler 1999). Once the mass-luminosity relation is 



obtained, the mass (and redshift)-flux relation trivially follows. 
5.2. SZ scaling relations 

The SZ effect is a scattering process appearing in the CMB 
spectrum as absorption at frequencies below ~ 218 GHz and as 
emission above. Nonetheless, a flux density can be formally as- 
sociated with the temperature distortion imprinted by the ther- 
mal SZ effect in the following way. 

We consider the Compton-y parameter observed in a given 
direction 8 of the sky and integrate it over an arbitrary solid 
angle £2 to obtain 



/ 

J a 



y(ff)d 1 6. 



(26) 



The temperature distortion over the patch of the sky covered 
by Q is proportional to Y times the typical frequency pattern 
of the thermal SZ effect, hence the monochromatic SZ flux per 
unit frequency received from the solid angle Q. can be defined 
as S v = j(v)Y, with 



j(y) = 2 



(kT y f 
{he) 2 



\m\ ■ 



(27) 



Here, T y is the CMB temperature and /(v) is the typical spec- 
tral signature of the thermal SZ effect, 



/(v) = 



(e x - l) 2 



+ 1 



1 



-4 



(28) 



where x = hv/kT y , and relativ istic corrections are ignored (see 
however Itoh & Nozawa| [2004 and references therein). 

In lSehgal et al.l (12007b . a prescription for linking the mass 
of a cluster to the Compton parameter integrated over the solid 
angle subtended by the virial sphere is proposed based on nu- 
merical simulations, according to 



*200(^200»Z) 



2.504 x 10" 



(D A (z)/lMpc) 2 \10 15 M, 



M 20 o 



1.876 



E(z) 2, \ (29) 



If we wish to convert the relation in Eq. d29l into a scaling 
law for the SZ equivalent flux density, the slope and redshift 
dependence will obviously remain unchanged. The normaliza- 
tion, however, must be converted into an SZ flux. We shall as- 
sume that the SZ monochromatic flux relates to I200, in the 
sense that the amount of SZ signal detected outside the virial 
radius is negligible. Keeping in mind that 



} (kT y ) 3 
' (he) 2 



2.701 x 10~ 18 J irT 2 = 2.701 x 10 n mJy, 



(30) 



the scaling law of Eq. d29l turns into 

/ M 20 o \ 

^ 2 \io 15 mJ 



(D A (z)/lMpc) z 

As stated in Sect. [4] for the SPT catalog construction we as- 
sumed v = v = 150 GHz, for which /(v ) = -3.833. The 
negative value indicates absorption, as is to be expected since 
v < 218 GHz. It follows 



S Va (M 2 oo,z) = 



2.592 x 10 8 mJy / M 2 oo 



1.876 



E(z) 2n . (32) 



(D A (z)/lMpc) 2 '\10 15 M, 
We use this scaling relation for the SPT cluster catalog. 

6. Catalog properties 

Figure [2] shows the minimum mass as a function of redshift 
for a cluster to enter each catalog, for all the dark-energy mod- 
els employed in this work. To compute that, we simply con- 
verted the limiting flux or integrated Compton-y parameter for 
the various surveys into a mass by means of the scaling re- 
lations described in Sect. [5] We notice that, as expected, the 
minimum mass for X-ray selected catalogs is a monotonically 
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Fig. 2. The minimum mass for clusters included in the catalogs 
constructed with the five different surveys investigated in this 
work, as labelled in the plot. Results for the seven different 
dark-energy models described in the text are shown, using the 
same color and line types as in Fig.Q] 

increasing function of redshift, while this is not the case for 
SZ catalogs, whose minimum mass slightly decreases at high 
redshift. This because the intrinsic redshift-independence of the 
SZ decrement causes the SZ flux density or integrated Compton 
y-parameter to scale as the inverse of the angular-diameter dis- 
tance squared, while the X-ray flux scales as the inverse of the 
luminosity distance squared, and the former tends to flatten at 
high redshift. An exception to this behavior occurs for Planck, 
whose limiting mass behaves in a more similar way to the X- 
ray catalogs because Planck has a very large beam that signifi- 
cantly smoothes the signal, especially when the angular size of 
the source is small. 

The differences between different cosmologies arise be- 
cause the scaling relations used in this work include distances 
and the expansion history of the universe, and hence depend 
on cosmology. The justification for this is that in models with 
dynamical dark-energy, and in early dark-energy models in par- 
ticular, structure formation begins at earlier times than in more 
standard models with w x = constant. As a consequence, clus- 
ters at a given redshift have more concentrated host dark-matter 
halos and more compact gas distributions, which enhances the 
SZ effect and X-ray emission. This agrees with the differences 
in the minimum mass for a given cosmology being more pro- 
nounced for X-ray catalogs, because X-ray emission is propor- 
tional to the square of the gas density, while the SZ effect scales 
only linearly with the density. 

Among the X-ray cluster catalogs, XCS has a systemat- 
ically higher minimum mass compared to eRosita, because 
of the higher minimum observed flux. Similarly, ACT has a 
higher minimum integrated Compton y-parameter than SPT, 



Fig. 3. The redshift distribution (all-sky equivalent) of clusters 
entering each of the five catalogs used in this work. Differences 
between different cosmologies are shown with the same color 
and line types as in previous figures, as labelled in the plot. 

and hence the minimum mass included in the catalog is sys- 
tematically higher. The minimum mass for the Planck catalog 
is always much larger than that of SPT and ACT, thus we ex- 
pect the number of clusters per unit area entering in its catalog 
to be relatively small. However, this is in some way compen- 
sated by the large area of the Planck survey. 

The redshift distribution of objects entering the various 
cluster catalogs for the seven different cosmologies used in this 
work is shown in Fig. [3] As expected, large differences oc- 
cur between different survey catalogs and different dark-energy 
models. The models EDE1 and EDE2 show very similar re- 
sults, and also the largest number of objects included in a 
given catalog. This is due to the enhanced c luster abundance in 
these models dFedeli & Bartelmann l2007ah . The model EDE3 
shows results that are very similar to the model with con- 
stant w x = -0.8, except for a moderate difference at low red- 
shift s, close to the peak o f the di stribution. This is agrees with 
Waizmann & Bartelmann! J2008). where it was shown that the 
cluster number counts predicted to be obtained with Planck in 
model EDE3 is almost identical to those for a standard ACDM 
model. The ACDM model always produces the smallest num- 
ber of objects, for all catalogs and at all redshifts, with the 
K08 model being intermediate between the ACDM and the 
w x = -0.8 models. 

As expected, the SZ surveys (except for Planck, which has 
a minimum mass behavior more similar to the X-ray catalogs) 
have a far wider redshift distribution than the X-ray catalogs. 
This is because the minimum mass slightly decreases at high 
redshift, as opposed to a monotonic increase (see Fig. |2j, and 
the sample may contain a high number of low-mass objects. 
For SPT we caution that the redshift distribution remains sig- 
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Fig. 4. Effective bias for the five cluster catalogs used in this 
work. Different colors and line styles refer to different cosmo- 
logical models, as labelled in the plot. 



nificantly above zero at the limiting redshift of our analysis, 
Z = 3. However, as we verified, the number of objects included 
in this catalog with z > 3 is just ~ 0.4% of the total, and hence 
negligibly contributes in the redshift integrals needed for com- 
puting the observed spatial and angular correlation functions. 
Those objects might be significant when binning the catalog 
in redshift, but again, the number of clusters with z > 3 is at 
most a few per cent of the number with z > 1.5, hence we as- 
sume that those objects can be safely neglected. The scaling 
relation between the SZ flux density and the mass of the host 
dark-matter halo is indeed highly uncertain for low masses and 
high redshifts, and hence we prefer to cut our sample at z = 3. 

7. Results 

7.1. Full catalogs 

In Fig. H] we show the effective bias computed by using all 
the clusters in the different samples investigated in this work. 
Results for the seven different cosmological models described 
in Sect. [2] are also shown. The effective bias for models EDE1 
and EDE2 is always significantly smaller than for the other 
models, and this holds true for all cluster catalogs considered 
here. This is obviously because forming massive objects is eas- 
ier in those models, and a much higher abundance of objects 
is present in the various catalogs than for the other cosmo- 
logical models. As a consequence, large galaxy clusters are 
less exceptional objects, and are less biased with respect to 
the underlying dark-matter density field. Models EDE3 and 
EDE4 (especially the former) are more similar to the standard 
ACDM case with respect to bias. For EDE3, this is consis- 
tent with the previous discussion and also with the findings of 



Waizmann & Bartelmannl Q2008), while for the EDE4 model 
this is unexpected, since the cluster abundance at a given red- 
shift (see Fig. |3]l is lower than, but quite similar to, those for 
EDE1 and EDE2. However, model EDE4 has an extremely low 
normalisation for the power spectrum of linear density fluctua- 
tions, which is expected to produce a higher 'monochromatic' 
bias. This is likely to play a major role in the computation of 
the effective bias. 

In line with the previous discussion, the ACDM model pro- 
duces the largest bias, while the K08 and w x = -0.8 mod- 
els give slightly smaller and very similar results to EDE3 and 
EDE4. In particular, the bias for the EDE1 model is up to a 
factor of ~ 2 smaller than that for ACDM at high redshift. It 
is important to note that, even though this is not clearly visible 
in Fig. |U the ratio of the effective bias in models EDE3 and 
EDE4 to that in model ACDM actually decreases quite steeply 
with redshift for z < 0.2, and increases again at higher redshift. 
Conversely, for other EDE models this ratio continously in- 
creases. This is due to the peculiar behavior of the dark-energy 
equation-of-state parameter in models EDE3 and EDE4, and 
will have important consequences on the effect of binning the 
cluster catalog in redshift, as discussed in the next subsection. 

In Fig. [5] we show the observed correlation function for the 
five surveys and seven dark-energy models considered here. 
As expected, the spatial correlation function decreases with 
increasing radius, and starts oscillating for large separations, 
r > 40 Mpc ft -1 . The correlation function in the standard 
ACDM model is practically indistinguishable from that in the 
model with a constant dark-energy equation-of-state parameter 
w x = -0.8 and also in the dynamical-dark energy model K08. 
In both cases, this is due to the similar behavior of the dark- 
energy equation-of-state parameter (see Fig. U}. In more detail, 
for the K08 case w x differs from the concordance value only at 
very low redshift, but approaches -1 at z > 3. In particular, the 
equation of state parameter is identical to the ACDM scenario 
during the linear stage of structure formation. 

On the other hand, models with early dark-energy differ sig- 
nificantly. Because of their lower effective bias, the observed 
correlation functions are also lower than those for more stan- 
dard models, by ~ 50 - 60% at small radii for eRosita and SPT, 
and slightly less for the other surveys. These observed differ- 
ences are caused by a combination of different effective bias, 
linear density-fluctuation correlation functions, and object red- 
shift distributions. 

The error bars in Fig. |5](shown only for the ACDM model 
for simplicity) were computed with the bootstrap method and 
imply that the difference between EDE and ACDM-like models 
would be detectable in the correlation functions observed with 
eRosita, SPT and (to a lesser extent) Planck. On the other hand, 
this difference would be completely lost in the noise for ACT 
and in particular for XCS. The large error bars visible in the 
XCS panel are due to the very small sky coverage of this survey, 
which is not adequately compensated by an increase in depth . 



This result qualitatively agrees with iMoscardini et al. I d2000h 



who showed that a deep survey has larger errors in the observed 
and angular correlation functions than a wide one. The same 
line of reasoning also applies to ACT, whose limiting integrated 
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Fig. 5. Predicted correlation function for the five cluster catalogs used in this work. Different colors and linestyles refer to different 
cosmological models, as labelled in the plot. Errorbars are computed with the bootstrap method, and refer only to the ACDM 
model for clarity. 



Compton->' parameter is too shallow to allow the collection of 
a significant signal. 

We note that in an attempt to reduce the size of the error 
bars, we increased the radial binning for XCS and ACT, so that 
the relative error, scaling oc 1/ y/dr, would decrease. However, 
even then the size of the error bars remains much larger than the 
differences in correlation amplitudes between the EDE models 
and the concordance cosmological scenario. We also note that 
in principle the same procedure could be applied to the rela- 
tive errors in the other catalogs, resulting in an even lower am- 
plitude of the relative errors. This is obviously unnecessary in 
this case, and in the following we also perform this operation 
only for the XCS and ACT catalogs. However, we keep in mind 
that when the error bars for the other catalogs are only slightly 
larger than the difference between cosmological models, en- 
larging the radial binning would probably allow a significant 
detection of deviations from the concordance ACDM model. 

We note that while Planck and XCS have a similar all-sky 
equivalent redshift distribution of objects, the sizes of the rela- 
tive errors are very different. This is because of the differences 



in sky coverage of about three orders of magnitude between 
the two different surveys, that enter quadratically in the com- 
putation of the error bars (see also the comments presented in 
Sect.[6l). 

A popular measure of the clustering amplitude from a given 
cluster catalog is the correlation length, defined as the spatial 
separation ro for which the correlation function equals unity, 
i.e., ^obst/o) = 1- Among the five catalogs analyzed here, the 
one with the largest correlation length is ACT, having ro — 16 
Mpc hT x for the ACDM cosmology, while the one with the 
smallest correlation length is eRosita, with ro - 10 Mpc hr x 
for the same model. As a consequence of the smaller correla- 
tion amplitude, the correlation length is also smaller in models 
with an early dark-energy contribution, by ~ 20% compared to 
the concordance cosmology. 

If the redshift information about clusters in one catalog is 
inaccessible or inadequate, i.e., if one has only projected infor- 
mation on the plane of the sky, then the only accessible clus- 
tering measure is the angular correlation function discussed in 
Sect. |3] Figure [6] shows the angular correlation functions. The 
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Fig. 6. Angular correlation function for the five cluster catalogs used in this work. Different colors and linestyles refer to different 
cosmological models, as labelled in the plot. Errorbars are computed with the bootstrap method, and refer only to the ACDM 
model for clarity. 



differences between different cosmological models and differ- 
ent surveys are enhanced, probably because the angular corre- 
lation functions are integrals over the spatial correlation func- 
tions along the line-of-sight. This results in a slight separation 
between the ACDM model and the K08 and w x = -0.8 models, 
which is absent from the observed three-dimensional correla- 
tion function. Accordingly, in the angular correlation function 
the ratio of the early-dark energy models to the ACDM-like 
models can be as high as a factor of ~ 3. Finally, differences be- 
tween individual EDE models are also enhanced, showing that 
the cosmology producing the lowest angular-correlation ampli- 
tude at all separations and for all catalogs is EDE4. 

As for the observed spatial correlation function, bootstrap 
error bars are shown for the ACDM case only for clarity. It 
is apparent that the difference between ACDM cosmology and 
EDE models is detected significantly by eRosita, while with 
SPT and Planck this is likely only for the most extreme model 
EDE4, at least with the radial binning we used in the error com- 
putation. If this radial binning is increased, the differences be- 
tween other EDE models and the concordance cosmology be- 



come detectable. Finally, as for the spatial correlation function, 
the error bars for XCS and ACT are too large to allow signifi- 
cant detection of any difference between cosmologies. In both 
cases, the survey is not deep enough to allow a noteworthy re- 
duction in Poisson noise. 

7.2. Redshift selected catalogs 

In the best-case scenario, knowledge of the redshift of objects 
entering the catalogs can be assumed. If this is the case, the 
galaxy-cluster sample can be sub divided into a limited num- 
ber of redshift bins, such that the number of observed (spatial 
and/or angular) pairs of objects is broadly the same in each 
bin. This has the effect of making the relative errors in the cor- 
relation functions approximately similar in each bin, and hence 
allows the coherent study of the redshift evolution in the cor- 
relation function. Measuring the redshift evolution in the clus- 
tering properties of galaxy clusters on the observational side 
allows us to determine more accurately the underlying cosmol- 
ogy and dark-energy evolution, since it provides an additional 



C. Fedeli et al.: Clustering with dynamical-dark energy 



11 




Fig. 7. Same as Figure[5]but with cluster catalogs restricted to z < 0.1 (left panel) and z > 0.1 (right panel). For the SPT catalog, 
the right panel actually shows the restriction to 0.1 < z < 0.3. Only three cosmological models are shown here for simplicity. 



constraint. As a byproduct of this analysis, we can also check 
whether suitable redshift binning can increase the ratios of the 
correlation amplitudes in the different cosmological models. 

We note that similar numbers of angular pairs in each red- 
shift bin should correspond to similar absolute number of ob- 
jects per bin, while this is not necessarily the case for spatial 
pairs. 

We then assume a perfect redshift knowledge of the ob- 
served sample. Tests showed that a redshift binning that en- 
sures a suitably high total number of clusters per bin is z < 0. 1, 
0.1 < z < 0.3, and z > 0.3 for all catalogs. This choice yields 
approximately equal signal-to-noise ratios in each bin for the 
angular correlation function. The same could be obtained for 
the observed spatial correlation function if each of the bins 
contained the same number of three-dimensional pairs. This is 
not true because the highest-redshift bin contains too few spa- 
tial pairs compared to the other two. The SPT catalog is the 
only exception due to its wide redshift distribution (see Fig. [3j. 
Thus, we shall use the same redshift binning for the spatial cor- 
relation function only for SPT and two bins, z < 0.1, z > 0.1, 
in all other cases. 

Figure [7] shows the spatial correlation function for all clus- 
ter catalogs studied here, for z < 0.1 in the left panel and 
z > 0.1 in the right. For the SPT catalog, the right panel 
shows results for 0.1 < z < 0.3, while results for z > 0.3 
are shown in Fig. [8] Here, we consider only three cosmolog- 
ical models, namely the standard ACDM and the two early- 
dark energy models EDE1 and EDE4. This choice was made 
because the EDE4 and ACDM models exhibit the largest dif- 
ferences in terms of spatial and angular correlation functions. 
Additionally, model EDE1 was chosen as being representative 



cr ease with increasing reds hift, in agreement with the findings 
of iMoscardini et al.l (120001) . The increment is significant for 



of early-quintessence cosmologies with a w x (z) evolution that 
completely differs from EDE4. 

The spatial correlation functions computed in all cat- 
alogs and for the three different cosmological models in- 
creasing red 
i et all d200( 

eRosita and Planck, and can become significant for SPT if the 
radial binning is increased, according to the previous discus- 
sion. For example, the correlation length for Planck increases 
from ro - 9 Mpc h~ x to ro ^ 17 Mpc hr l in going from z < 0.1 
to z > 0.1. As for the SPT catalog, for which we added a high- 
redshift bin, the correlation length increases from ro - 7 Mpc 
/T 1 forz < 0.1 up to r - 14 Mpc h~ l forz > 0.3. 

The ratio of the correlation amplitudes in different cosmo- 
logical models generally increases when clusters are restricted 
to z < 0.1, and the ratio of the ACDM and EDE4 models in 
terms of the spatial correlation function increases significantly. 
For example, in the eRosita sample the ratio of the correlation 
functions increases from a factor of ~ 60% to ~ 2, and similar 
increments are seen in the other catalogs. On the other hand, the 
relative difference between EDE1 and ACDM is practically un- 
changed compared to the full samples. For high-redshift sam- 
ples, the ratio of the correlation amplitudes in ACDM and 
EDE4 models is still larger than for the complete samples, but 
smaller than for low-redshift samples. The differences between 
EDE1 and the concordance ACDM models are again very sim- 
ilar to those between the low-redshift and the full samples. 

The different behavior of EDE1 and EDE4 can be attributed 
to the different trend with redshift of the effective bias dis- 
cussed above. Due to the behavior of the redshift distributions 
shown in Fig. [3] the subsamples with z < 0.1 and z > 0.1 are 
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Fig. 8. As Fig. [5] but only for the SPT catalog restricted to z > 
0.3. Only three cosmologies are shown here for simplicity. 

virtually dominated by objects at z ~ 0.1 and z ~ 0.2 respec- 
tively (because of a combination of a sharp decline in the mass 
function with redshift and volume effects). The ratio of the cor- 
relation functions of the underlying density fluctuations in the 
EDE and ACDM models are virtually unchanged when going 
from z ~ 0. 1 to z ~ 0.2, and the same is true for the ratio of the 
effective bias in ACDM and EDE1 model, although in the lat- 
ter model it increases steadily, but slowly with redshift. Instead, 
the ratio of the biases in ACDM and EDE4 decreases steeply 
between z ~ and z ~ 0.2; The difference between the EDE4 
and ACDM model is therefore expected to be more enhanced 
in low -redshift than high-redshift catalogs, while the difference 
with EDE1 remain about the same. 

The relative errors are still very small for eRosita in all 
redshift-selected catalogs, thus the signal-to-noise ratio for dif- 
ferences in the correlation amplitudes between ACDM and 
EDE models is optimal for the z < 0.1 catalog. The same is 
also true for Planck, because the error bars tend to increase in 
the z > 0.1 sample. For SPT, the only redshift bin in which a 
difference between the ACDM and the EDE 1 or EDE4 models 
can be reliably detected is the high-redshift bin, z > 0.3. For 
this bin, we obtain essentially the same results as for the full 
sample. The situation remains practically unchanged for XCS 
and ACT, where the relative errors are still very large and do 
not allow any significant detection of the differences between 
concordance and more exotic models. 

Figures [9] and [10] show the angular correlation functions 
for cluster catalogs binned in redshift according to the scheme 
Z < 0.1, 0.1 < z < 0.3 andz > 0.3, as discussed above. For con- 
venience, we show results for all the catalogs considered here 
and for the cosmological models EDE1, EDE4 and ACDM. In 
contrast to the spatial correlation functions, the angular correla- 



tion functions decrease in amplitude with increasing redshift, a 
trend that is significant for the eRosita and SPT samples, while 
the large relative errors for the Planck catalog probably allow 
a detection of this decrement only at large angular separations 
and if the radial binning for the computation of error bars is 
enlarged. 

The ratio of the correlation amplitudes in the ACDM model 
to those in the EDE models is slightly lower in the first two 
(low-redshift) bins than for the full samples. Since the relative 
error bars also increase in those bins, the signal-to-noise ratio 
for differences in the angular correlation functions is also lower. 
In contrast, the ratio increases slightly in the highest-redshift 
bin, z > 0.3, but there the relative errors are also larger, espe- 
cially for Planck. As a consequence, differences between the 
concordance model and both EDE models can be significantly 
detected with the eRosita and, at least at large angular separa- 
tions, with the SPT full sample or selected in the intermediate 
redshift bin, 0.1 < z < 0.3. The situation remains again un- 
changed for XCS and ACT, with relative errors being too large 
to reach any significant conclusion. 

According to the discussion above, distinguishing between 
different cosmologies by means of the angular correlation func- 
tion does not take any particular advantage of redshift selection. 
It is legitimate to ask if the situation is different when selection 
is performed based on the observable used for cluster survey, 
instead of pure redshift. For instance, binning X-ray selected 
clusters according to their flux would correspond to a simulta- 
neous (and highly non-linear) binning with mass and redshift. 
To check this issue, we performed clustering analysis of the 
SPT cluster sample after binning it according to flux density as 
5 mJy < S Va < 10 mJy, 10 mJy < S Vo < 25 mJy and S Vo > 25 
mJy. This choice ensures a suitably high number of spatial and 
angular pairs in each bin. However, the outcome is that the error 
bars are always larger than the differences between the various 
models for all bins. In other words, this kind of binning does 
not improve the detection of deviations from the concordance 
ACDM cosmological model with respect to the angular corre- 
lation function, and hence we decided not to show these results 
here. 

8. Summary and discussion 

We have studied the clustering properties of galaxy clus- 
ters in various cosmological models with dynamical and non- 
dynamical evolution in the dark-energy component. In addi- 
tion to the concordance ACDM cosmology, we addressed a 
model with a constant equation-of-state parameter for dark en- 
ergy, w x = -0.8, a dynamical-d ark energy model wit h the 
w x (z) parametrization proposed by iKomatsu et alj (|2008), and 
four models with a non-negligible amount of early dark en- 
ergy. Cosmological models with dynamical dark energy are 
generically expected to form structures earlier, affecting both 
the abundance of massive galaxy clusters and their spatial dis- 
tribution and clustering properties, this latter property having 
been the issue here. 

To predict forthcoming observations, we computed the ef- 
fective bias as a function of redshift and observed spatial and 
angular correlation functions as a function of (physical or ap- 
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Fig. 9. As Fig. [6] but with cluster catalogs restricted to z < 0.1 (left panel) and 0.1 < z < 0.3 (right panel). Only three cosmological 
models are shown here for simplicity. 



parent) separation that are expected to be measured in cluster 
catalogs produced by planned blind surveys both in the X-ray 
and in the sub-mm regimes by means of the thermal SZ ef- 
fect. The X-ray surveys that we considered in this work are 
the eRosita wide survey, which is described in detail in the re- 
lated mission definition document, and the XMM cluster sur- 
vey, based on existing pointings of the XMM satellite. For 
the SZ surveys, we focused on the South Pole Telescope, the 
Atacama Cosmology Telescope, and the Planck all-sky survey 
of the Cosmic Microwave Background. 

For computing the clustering properties of objects con- 
tained in each catalog as a function of cosmology, we employed 
a well-established formalism that takes past-light cone and se- 
lection effects on the cluster sample into account. To link the 
limiting flux of each survey to the minimum mass that enters 
the respective catalog at a given redshift, we adopted realis- 
tic scaling relations, based both on observations and numerical 
simulations, between the mass of a cluster and its X-ray lumi- 
nosity or SZ flux density, or its integrated Compton-y parame- 
ter. It turns out that the minimum mass entering a cluster cata- 
log depends not only on the instrument considered, but also on 
the cosmological model adopted. This is so because the scal- 
ing relations mentioned above depend on the underlying world 
model through the expansion rate and cosmological distances. 

As one could naively expect, the number of objects entering 
a given catalog at a given redshift depends heavily on the cos- 
mology, the highest cluster abundances being present in models 
with EDE. This is caused by the well known higher mass func- 
tion displayed by these models and the lower minimum mass 
entering the catalogs. The SPT catalog is the most extended in 
redshift due to the low limiting SZ flux density (5 mJy), while 



distributions for XCS and Planck are the most limited in red- 
shift, because of the quite shallow limiting X-ray flux of the 
former and large beam of the latter, which tends to dilute the 
signal. 

For all catalogs, the first two models with early 
quintessence, EDE1 and EDE2, display the smallest effective 
bias at all redshifts. This is due to the high abundance of struc- 
tures present in these models. On the other hand, the concor- 
dance ACDM model always displays the highest effective bias, 
for all catalogs and at all redshifts, while the other models lie 
somewhere in-between the two. Because of the width of their 
extent in redshift, SPT and ACT catalogs show the flattest trend 
for the effective bias, especially SPT, for which b e ff(z) is at most 
~ 10 at z = 3. Other catalogs, particularly those based upon 
Planck and XCS have a much steeper trend, the effective bias 
reaching considerably higher values already at z ~ 2. 

Concerning the spatial correlation function, all the EDE 
models almost coincide on all scales, having a correlation am- 
plitude smaller than that for the models ACDM, K08 and con- 
stant w x = -0.8, which are also very similar to each other. 
The largest difference is displayed by the catalogs produced 
with eRosita and SPT. In the former, the relative errors are also 
smallest, thus making eRosita the most promising instrument 
for the detection of EDE through the spatial correlation func- 
tion, if the full cluster catalog is to be used. Detection might 
also be possible for SPT and Planck at intermediate scales 
r ~ ro, while it is completely out of question for XCS and 
ACT, whose error bars are too large. This is due to the shallow 
limiting fluxes for both of them, which is not adequately com- 
pensated for by the area covered (as is the case for Planck). The 
situation is similar for the angular correlation function. There 
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the smallest error bars are produced by eRosita as well, while 
the EDE4 model produces the smallest correlation amplitude 
for all the catalogs. 

To explore the redshift evolution in the spatial and angular 
correlation functions, we also performed different cuts in red- 
shift of the various cluster catalogs, in a way that the number 
of pairs of objects in each bin would remain approximately the 
same. For the number of spatial pairs, we adopted the double 
binning z < 0.1 and z > 0.1, with the exception of SPT which 
allows the inclusion of a third bin, z > 0.3, with the modifica- 
tion of the second bin to 0.1 < z < 0.3. This same binning was 
also employed for the number of angular pairs. 

The result of this analysis is that the spatial correlation 
function increases with increasing redshift, the correlation 
length for eRosita growing by ~ 80% between the low- and 
high-redshift bins. The relative errors show that the increment 
in the correlation function is significant for eRosita, SPT and 
Planck, for all cosmological models considered in this work. 
Also, comparing the amplitude of the ratio between spatial cor- 
relation functions in EDE models and ACDM-like cosmolo- 
gies, with the size of the error bars indicates that it is better to 
focus on the low-redshift cluster subsample in order to max- 
imize significant differences between the concordance model 
and the EDE models in the eRosita and Planck catalogs. The 
high-redshift bin, z > 0.3, is better for the SPT catalog, giv- 
ing results compatible with those from the full catalog. As for 
the angular correlation function, it tends to decrease with in- 
creasing redshift, a trend that is significant for all models and 
catalogs except the usual XCS and ACT. In the absence of er- 
rors, in the high-redshift bin there would be the highest chance 
of distinguishing a ACDM model from models with an early 
quintessence contribution, since there the ratios of correlation 
amplitudes are at their highest. However, the errorbars are also 
large there, so that eRosita is the only survey expected to permit 
significant detection of deviations from the concordance model 
at all redshifts. In general, if one is interested in optimizing the 
differences between angular correlation functions measured in 
different models, it does not pay to subdivide catalogs accord- 
ing to redshift. 

We also found that the same is true when cluster catalogs 
are binned according to the observable used to define them, an 
approach more directly motivated from the observational point 
of view. Specifically, binning the SPT catalog according to the 
flux density always produces ratios of angular correlation func- 
tions in different cosmologies that are similar to those for the 
full catalog, while the errorbars are always larger. 

Before concluding, two notes of caution are in order. First 
of all, some discussion arose in the literature on whether 



the se mi-analytic calculations performed in Bart elmann et al 



(2006) on the spherical collapse model in cosmologies with 
dynamical-dark energy are indeed correct, and h ow accurately 



they a re reproduce d by jV-body simulations (iFrancis et al 



2008afllGrossi & Springe3E )08). While the discussion is not 
yet settl ed. | Sadeh et all d2007l) p erformed the same calcula- 
tions as iBartelmann et alj d2006l) using a different approach, 
and found resu lts that are perfectly cons istent with the lat- 
ter. In addition, Schafer & Kovama (2008) used the same ap- 
proach as Bart elmann et alJ 12006) to successfully evaluate the 
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Fig. 10. As Fig. [6] but only for catalogs restricted to z > 0.3. 
Only three cosmologies are shown here for simplicity. 

spherical collapse behavior in a modified gravity scenario. 
Hence, we a re at least reasonably con fident that the approach 
followed by Bartelmann et al. ( 20061) . also employed in this 
work, is fundamentally correct. Moreover, numerical simula- 
tions using the correct early-time behavior of the growth fac- 
tor for s caling the initial conditio ns y ield results o t her tha n 
those of lGrossi& Springell d2008l) and IFrancis et all d2008ah . 
which tend towards the expectation from the analytic work of 
Bartelmann et al. ( 2006h . Although precise direct integrations 
of the spherical collapse equations are difficult, further work 
avoiding approximations has so far confirmed our earlier re- 
sults (Pace et al., in preparation). Even though definitive con- 
clusions are not reached yet, these facts seem to justify our 
choice. 

Secondly, estimates of the number counts of galaxy clusters 
detected by Planck seem to fal l subst antially below the esti- 
mates of I Schafer & Bartelma nn (2007), and the redshift distri- 



bution is apparently shallower than that represent ed in Fi g. [3] 
(Plan ck SZ Challenge (in preparation), see also iLeach et alj 
2008). While definitive new limits for Planck are not yet avail- 
able, this part of the results should be read with caution. In 
particular, a decrease in the number of objects in the Planck 
catalog would produce an increase in the relative errors, which 
scale as the inverse of the square root of the number of pairs of 
objects. 

The present work shows that, while mild modifications to 
the redshift evolution in the equation-of-state parameter for 
dark energy w x {z) have a negligible impact on the clustering 
properties of galaxy clusters, more exotic models such as early 
dark-energy cosmologies can change the effective bias of col- 
lapsed objects significantly, and thus also the spatial and angu- 
lar correlation amplitudes of galaxy clusters. We have shown 
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that at least some of the forthcoming blind surveys both in the 
X-ray and sub-mm regimes will be able to distinguish signif- 
icantly between these models and more generally place con- 
straints on the time evolution in the dark-energy density. 

In general, we expect object clustering to be less effective 
than e.g., direct abundance data in determining the cosmolog- 
ical model. For instance, in the eRosita catalog the number 
of clusters at z > 1 increases by about one order of magni- 
tude between the ACDM and the EDE4 models. This variation 
is much larger than the corresponding variation in the corre- 
lation functions, while the size of the relative errors, assum- 
ing Poisson statistics, is comparable. Object counting at high- 
redshift is also expected to be capable of distinguishing cosmo- 
logical models with a gentle variation in dark-energy density 
from standard cosmology. The number of high-z clusters in the 
eRosita catalog is higher by a factor ~ 2.5 in the model with 
constant w x = -0.8 than for ACDM. 

Nevertheless, the results of this work show that clustering 
of massive clusters by itself remains a fundamental channel to 
unravel the effect of the expansion history of the Universe on 
the process of structure formation, although employing this in- 
formation in addition to simple object number counts is cer- 
tainly an interesting issue to be explored. 
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